arXiv:1509.07697v2 [nucl-th] 17 May 2016 


Matching effective few-nucleon theories to QCD 

J. Kirscher 

Racah Institute of Physics, The Hebrew University, Jerusalem 91904, Israel 

The emergence of complex macroscopic phenomena from a small set of parameters and microscopic 
concepts demonstrates the power and beauty of physical theories. A theory which relates the wealth 
of data and peculiarities found in nuclei to the small number of parameters and symmetries of 
quantum chromodynamics is by that standard of exceptional beauty. 

Decade-long research on computational physics and on effective field theories facilitate the as¬ 
sessment of the presumption that quark masses and strong and electromagnetic coupling constants 
suffice to parameterize the nuclear chart. By presenting the current status of that enterprise, this 
article touches the methodology of predicting nuclei by simulating the constituting quarks and glu¬ 
ons and the development of effective field theories as appropriate representations of the fundamental 
theory. 

While the nuclear spectra and electromagnetic responses analyzed computationally so far with 
lattice QCD are in close resemblance to those which intrigued experimentalists a century ago, 
they also test the theoretical understanding which was unavailable to guide the nuclear pioneers 
but developed since then. This understanding is shown to be deficient in terms of correlations 
amongst nuclear observables and their sensitivity to fundamental parameters. By reviewing the 
transition from one effective field theory to another, from QCD to pionful chiral theories to pionless 
and eventually to cluster theories, we identify some of those deficiencies and conceptual problems 
awaiting a solution before QCD can be identified as the high-energy theory from which the nuclear 
landscape emerges. 
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FIG. 1. (Color online) Range of applicability of nuclear effective theories (adapted from Ref. [4] ) in a space spanned by the 
pion mass m,r, a low-momentum scale Q, and the nucleon number A. The physical world is measured at ~ 140 MeV on the 
red sheet. LQCD assesses the few-nucleon systems of the deuteron d, the triton t, and the a-particle on the dashed enclosed 
sheets at 510 and 806 MeV. EFT(/) (light green) was applied to ^4 < 5 nuclei up to large pion masses. XEFT (yellow area) 
breaks down at some m^, but applies at larger momenta than its perturbative version and EFT(^). Cluster EFTs (cyan) have 
only been used at physical m,r. 


I. INTRODUCTION 


Lattice QCD as an experiment with light baryons — It is a striking analogy that the calculation of binding energies 
of the lightest nuclei as compound objects, whose constituent dynamics are dictated by quantum chromodynamics 
(QCD), proceeds along similar paths as the early experiments in the real world. The relative determination of the 
deuteron mass to Helium mas^, with its preceding preparation of a sufficiently pure beam of the respective nuclei, 
resembles remarkably the lattice methodology. In the latter, mass ratios are also more accessible because of a reduced 
statistical uncertainty. Furthermore, devising interpolating operators with a large initial overlap with the baryons of 
interest are the numerical analog of preparing a clean beam. The exploration of the nuclear landscape included in 
parallel to the study of mass spectra the response of nuclei to electromagnetic field^. Again, these investigations 
find their resemblance in contemporary lattice calculations of the magnetic structure of light nuclei. The fact that 
the lattice simulates nuclei at artificially large pion masses yields more than just data points for an extrapolation to 
the physical pion mass. It allows for the investigation of phenomena emergent from a theory as sophisticated as the 
one we deem appropriate for the strong force in nature. The theory at unphysical pion masses shares symmetries and 
degrees of freedom and thereby can identify features that are peculiar to QCD as the theory which breaks the flavor 
symmetry with some specific quark-mass values or can be universally attributed to the three SU(3) color and flavor 
symmetric quarks. 

Tools that have evolved over the century and the questions which can be investigated with the lattice technology 
are plenty. The basic motivation, to build a theory on the available data to predict what presumably will be measured 
with more effort in nuclei with a total number of neutrons and protons A > 4, using different external probes, or 
looking at the systems at different scales, remains, however, invariant. Today, we start at the same point like physicists 
when facing the first measurements on light nuclei in the physical world. 

Contemporary theoretical nuclear reality — With the increasing resolution of experimental apparatus and the 
ensuing discovery of the substructure of nucleons the interaction amongst them lost its status of a fundamental theory 
and became to be thought of as a manifestation of QCD at low energy. The current understanding of nuclei comprises 
a set of effective theories, each devised to describe them only up to a certain nucleon number A or momentum scale 
Q, at which the neutrons and protons can be excited and/or deformed but their internal structure is not probed. 
The pioneering theories and models for the nuclear interaction, e.g. , meson-exchange (Yukawa) theories and quark 
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models, can be recoverecQ from or found to be included in these effective theories. 

Specifically, we understand the nuclear interaction to be given by chiral effective field theory (XEFT) which: (i) 
treats the pion-nucleon coupling perturbatively for very low momenta, (ii) becomes non-perturbative for higher 
energies, (Hi) can be matched to the pionless effective field theory (EFT(/)) which is applicable in the baryon sector 
for momenta much smaller than the nuclear scale set by typical binding momentum in nuclei of about 100 MeV/c. 
The latter serves itself as basis for cluster theories which introduce composite nuclei, like a particles or deuterons, for 
the description of larger systems. This understanding is at present incomplete as a reconciliation of differing power¬ 
counting schemes which define XEFT is yet to be found. A fully renormalization-group invariant version of XEFT, for 
example, has not yet been formulated. However, given a reliable counting scheme and thereby uncertainty estimates 
will enable us to relate the five relevant QCD and QED parameters (three light-quark masses and the electromagnetic 
and strong coupling constants) via the matching of QCD and EFT amplitudes to observables, which are poorly known 
experimentally and would thus be an ultimate assessment of whether or not the nuclear theory has been found. 

The crucial calculation of multi-hadron amplitudes with QCD became reality, recently, via the numerical solution 
of the QCD path integral in discrete space-time (lattice QCD or LQCD): unquenched, high-statistics LQCD mea¬ 
surements are available for the nuclear spectrum up to A = 4 using (i) Nj = 2 + 1, i.e., three-flavor QCD with the 
physical strange-quark mass and equal up, down-quark masses m„ = rud corresponding to ~ 300 MeV in Ref. [5] , 
450 MeV in Ref. [6] (A < 2), and 510 MeV in Ref. [7] , and (ii) Nj = 3, i.e., exact SU(3)-flavor QCD because of 
degenerate light-quark masses yielding m,r 806 MeV in Ref. [8] . While measuring at unphysically high pion masses 
increases the signal-to-noise ratio in lattice simulations and allows for meaningful results, it inhibits a matching of 
those amplitudes to a chiral effective theory which relies on pion masses < 500 MeV. Before LQCD can probe the 
nuclear spectrum in this pion-mass range, its amplitudes can be matchecplin the region where its area of applicability 
overlaps with the pionless theory. 

The current status of this program, to relate QCD parameters through LQCD, XEFT, and contact theories to few- 
nucleon observables, is reviewed in this article. A graphical outline is given in Fig. It shows areas of applicability 
of the various theories in the physical parameter space spanned by the pion mass, a momentum scale Q, and the 
nucleon number A. In this figure, QCD is shown to probe the entire space and is expected to coincide with experiment 
on the physical (red) sheet. Without QCD solutions for nuclear observables (dashed enclosed, transparent sheets) 
at physical m^n-, we begin the discussion at the intersection between experiment and XEFT (yellow sheet) in Sec, [h} 
before showing how the overlap between XEFT and EFT(^) was utilized to reach larger nuclei (green squares indicate 



include comments on the further investigation of the A axis by matching contact to cluster theories (cyan) in the 
outlook. 


II. LATTICE QCD ^ XEFT 

To predict LQCD observables at physical a theory is needed which is consistent with QCD at physical and 
some larger WTr, where it can be matched to available LQCD data. Furthermore, the power counting of that theory 
has to be understood for all values of the pion mass between the matching point and the physical point. The extension 
of chiral perturbation theory to multi-baryon systems, hereafter referred to as XEFT (concepts in Refs, [ini [H] , and 
reviews in Refs. [mini) is the only ansatz which has been used for such an extrapolation. The power counting 
of such an extensiorj^ however, is still under development (consider Refs. [TSIED] next to a solution suggested in 
Refs. [mill]). Given the development of such an effective field theory for few-nucleon systems and the availability 
of LQCD data at a where chiral perturbation theory (ChPT) converge^ the algorithm (original formulation in 
Ref. [H] ) reviewed below can be used to predict physical observables from LQCD data at unphysical 

The chiral Lagrangian depends explicitly and implicitly on the pion mass. The implicit dependence resides in the 
pion decay constant /^, the axial coupling constant gA, and the nucleon mass mpj. Using ChPT, this dependency can 
be made explicit in a power expansion, e.g. , for the decay constantP^lwith a renormalization scale g: 

Mm.) = Mm. = 0) (^1 + In ^ + 0(1)) + 0(m^)) . (1) 


^ See Ref. [3] as an example for a derivation of a quark model from the heavy-baryon effective field theory of QCD. 

^ Ref. |14| demonstrates the complexity of finding inconsistencies in a perturbative treatment of the pion exchange; Ref. |15| reviews also 
problems of other counting schemes. 

^ Ref. |23| hints to a breakdown at m-n- < 500 MeV of the chiral expansion for /tt. 
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The XEFT Lagrangian contains in addition to single-nucleon, single-meson, and meson-nucleon coupling parameters 
pion-mass independent and dependent contact operator^ 


>Cfour-fermi = Cs (N^N)^ + Ct {N^ 

+Dsi (N'^MN)'' {N'^N)+Dti {N^Ma,N)^ (N'^cr.N) 

+Ds2 (N^N)^ (N^N) TtM+ Dt 2 {N^a,N) TrM 


( 2 ) 


The associated low-energy constants (LECs) Cs,t and Di are pion-mass independenlQ They can be determined at 
any within the convergence radius of ChPT. From the available data at physical TOtt, Di cannot be determined 
independently from Ci. One way to resolve this ambiguity are additional scattering experiments of pions on deuterons 
or heavier targets. Another utilizes two-nucleon observables at different m^r as the necessary constrains to separate 
the explicit m,r dependence from the LECs as shown in Eq. pin One LQCD data set for a pion mass which might 
be in the range of applicability of XEFT was calculated in Ref. at TOtt = 450 MeV. 

Replacing LQCD data at another with experimental data to pin down the LECs, was demonstrated in Ref. m- 
The authors employed the Beane-Kaplan-Vuorinen (BKV) power-counting schem^^, a variant of the Kaplan-Savage- 
Wise (KSW) schem^^ for XEFT to relate the two-nucleon amplitude to the pion mass. The BVK scheme introduces 
in contrast to KSW an additional fictitious heavy meson of mass A. In effect, this particle cancels the r~^ singularity of 
the pion exchange in the triplet channel where KSW does not converge. The mechanism is similar to the Pauli-Villars 
regulation in QED. The scattering lengths at physical and at ~ 350 MeV of Ref. [SD] were used to renormalize 
the NN counter terms unambiguously. The pion-mass dependence of the effective range can be inferred from data at 
a single ttItt because up to NLO, operators proportional to the quark mass do not contribut^^. In addition to the 
ensuing dependence of the deuteron binding energy (as a pole of the amplitude Eq. ([^), Ref. [27] derived the 
charge radius, the magnetic moment and polarizability, and the photodisintegration of the deuteron as a function of 
the pion mass. We compare the binding energies in the singlet and triplet channels of Ref. m in Figj^ (thin red 
line) with the results obtained with a chiral potentiaPD. The latter neither binds the singlet state at a higher nij^, 
nor does it indicate limm^_,,o = 0. 

Before Ref. [27] , and without lattice data at sufficiently low pion masses, XEFT was used, also with a modified KSW 
counting which uses assumptions about details of the short-range structure of the interaction, for an extrapolation in 
in Ref. m- Like BVK, this Beane-Bedaque-Savage-van-Kolck (BBSvK) scheme is identical to KSW in the NN 
singlet channel but as an expansion about the chiral limit = 0 iterates the pion exchange in the coupled triplet 
channel (Weinberg counting). In Ref. [T^ , BBSvK is used to investigate the deuteron and ^S'o , neutron-proton 
(np) scattering properties for pion masses from 0 up to 200 MeV. The two-nucleon amplitude was matched at physical 
to the triplet np scattering length {^Unp) and effective range (ra) to determine the relevant contact terms in 
Eq. §. With a leading order chiral expansion for /^, g^, and tUn, the deuteron was found to become unbouncQat 
^ 100 MeV with a corresponding divergent ^a„p. 

The extension of the method to the np-singlet channeP^ demonstrated the relevance of the above ambiguity in the 
contact terms. One choice of Ci and Di left the di-neutron inn) unbound for all pion masses in [0,300 MeV] while 
another yielded nn bound in the interval between 180 and 240 MeV. The problem was then addressed in Ref. [33] by 
placing constrains on the contact terms with naive dimensional analysis. A relevant conclusion for the extrapolation 
from ttItt > 140 MeV down was the consistency of a bound nn, z.e., singlet two-nucleon state, with data at the physical 
pion mass and coupling strengths of natural size. 

The effect a different power counting has on the rn^ dependence of the two-nucleon system is seen when com¬ 
paring the results of Refs. [ 13131133 ] , to Ref. [ 33 ] , and to Ref. [ 33 ] . Considering also two-pion-exchanges in a 
modifiecP^ Weinberg counting. Ref. |34] found the deuteron in the chiral limit more strongly bound than in nature. 
While this is in contrast to the unbound deuteron of Ref. [15] , the absence of bound states in the singlet NN channel 
for TOtt G (0, 200 MeV) is a common result of both schemes, pions in Ref. [ 33 ] finds the singlet and deuteron state un¬ 
bound in the chiral limit. This work employs composite two-nucleon spin-singlet and triplet quasi-particles (dibaryons, 
see Sec. V B) coupled to single nucleons and propagates the pion-mass dependence to a pionless theory through the 
analytical matching of amplitudes. These differences should concern convergence rates and precision of the EFT but 
must eventually lead to the same observables. This issue is unresolved. As a feature, ref. [ 33 ] demonstrates a reversed 
extrapolation from the physical pion mass to larger values. The LQCD results which are available (see Sec. IV B) 


^ with M ~ ^diagirriu, m^) g = 1 -I- A— t ■ tt — + OItt®), Pauli isospin matrices t, pion isovector tt, and 

nucleon iso-doublet N = 

^ The pions are dynamical in XEFT. In EFT(]t), in contrast, the contact LECs do depend on mTr- 

^ The method of resonance saturation as developed in Ref. m and applied, e.g. , to the quark-mass sensitivity of a collection of nuclear 
observables relevant for big-bang nuclear synthesi^^, approaches the problem by relating XEFT contact LECs to phenomenological 
boson-exchange models. As it is unknown how the uncertainties of that method can be quantified, we abstain from further elaboration. 

^ The value was also founciSl to be quite sensitive to the regularization of the contact interactions. 
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FIG. 2. (Color online) Pion-mass dependence of the nn-singlet (top left) and deuteron (top right) binding energy, the neutron- 
deuteron quartet-channel scattering length (bottom right), and the neutron-deuteron(nn-singlet) (blue (red), bottom left) 
doublet scattering length. Results for m^< 220 MeV (left y scale) represent predictions with EFT(^) matched to XEFT (bound 


states are red solid line^^and red circleJ^, unbound intervals are gray lineJ^) and for m.,r> 230 MeV (right y scale) lattice 
data for B 2 and predictions of EFT(^) matched to LQCD data (see sections |IV B| and 
uncertainty and the vertical black line marks m.n- ~ 198 MeV at which —1 00 . A' 


V B I. Box heights correspond to lattice 
1 X scales are nonlinear. 


there provide support that nature as represented by experimental data at physical m,n. “flows” into SU(3)-large-TOq 
QCD as explored via the lattice method. 

Regardless of the counting scheme, a critical pion mass close to the physical value was found where the two- 
body scattering lengths diverge. The associated scale invariance of the two-nucleon system w.r.t. any coordinate 
transformation r —A r for A G K+ indicates an infrared fixed point of the renormalization-group (RG) flow of the 
QCD coupling constanl[^ How this infrared QCD fixed point is expressed in the RG flow of XEFT is unknown 
and presents a problem intimately related to the development of a consistent power counting. In contrast, the flow is 
known for an effective theory that considers solely contact interactions amongst nuclei (see Sec. |III| below). Although 
not manifest in its coupling constants, the peculiar limit-cycle trajectory of the leading three-body LEC, emergent 
from the fixed point, should find its signature in the AEFT spectrum. Part of this signature, namely remnants of an 
Efimov spectrunJ^, was found^when the authors employed an interaction derived in Ref. [34j in the three-nucleon 
system. The characteristic asymptotic ratio of binding energies at the deuteron-neutron accumulation point of ^ 515 
was approximated, while the ground state, be., the triton was stable = 4 ± 2 MeV) in the considered pion-mass 
interval between 190 and 210 MeV. The XEFT predictions in Fig. |^were obtained at a fixed cutoff and thus do not 
allow an assessment of the uncertainty like the displayed LQCD and physical experiments. 

In essence, all three approache^^^^^ I® employ XEFT to define a Lagrangian with implicit and explicit dependence 
on rTT-Tr. They differ in the power counting and thereby the dependence of observables on the parameters of £, be., m^. 
Pion-deuteron scattering experiments in the physical or lattice world separating the Di counter terms from the Ci’s, 
like (7V^7V)^7r7r, or LQCD amplitudes at different 111 ^^ will allow for a unique determination of the LECs, 

and thereby aid the development of the above mentioned consistent power counting. 


This nontrivial inference assumes a valid approximation of QCD by XEFT. Only then does a statement based on a XEFT calculation 
about a QCD parameter make sense. 
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To that end, the EFT analysis in Ref. [5] demonstrated a match of XEFT to LQCD instead of experimental data 
similar to Ref. [27] . In the KSW^and BBSvK?^ counting schemes, this work considered the EFT convergence in 
the ^S'o and ^Si scattering lengths and effective ranges up to NNLO with LQCD binding energies and phase shifts 
at TOtt = 450 MeV as input. The convergence rate was found small in all 4 observables and different from the ones 
extracted from the LQCD phase shifts via an effective range formula. The latter discrepancy was explainecPby fitting 
to data beyond the natural scale of the theory. 

The results at 450 MeV pion mass are entering a region where they can be analyzed with low-energy theorems as 
developed in Ref. [38j . These theorems provide expansions of the modified effective range parameters in 
This assumes a long-range interaction provided by XEFT, hence as the light scale, and a short-range potential 
of range m~^. The uncertainty due to an unknown m^n- dependence on the short-range part of the interaction was 
founcP^to increase significantly with niTr- Part of this uncertainty is, again, the isolation of the explicit pion-mass 
dependence in the four-fermi contact LEC^ The method applied to < 400 MeV yields results for a/r which 
are consistent within error bars with the m^n- ^ 450 MeV data (see Table |ll] and compare Fig. with Figs. 6 and 7 
in Ref. [38] ) regardless of being extrapolated from the physical m^n- or constrained by LQCD deuteron data directly 
at rriTr = 430 MeV. Even a naive continuation of the derived ttItt dependence of the ratio a/r is “in line” with the 
peculiar value of 2 LQCD finds at 806 MeV pion mass (see discussion around Eq. Q)- 

The above attempts employ effective field theories to parameterize two- and three-nucleon amplitudes, i.e., the 
scattering lengths in singlet and triplet S-wave channels and the triton binding energy with the pion mass. 

The pion mass and the theoretical uncertainty of XEFT can then be passed on to another effective theory which 
is more practical for larger nuclear systems. This matching between a pionful and pionless theory cited above in 
Refs. [57| |3S| , leads to the next section where we summarize work which already employed the idea to propagate 
the pion-mass dependence from two- into few-nucleon amplitudes with the caveat of an unknown uncertainty in those 
amplitudes and therefore itIt^ dependence of the underlying theory which is then no longer QCD but XEFT. 


III. XEFT ^ EFT(/) 


Above, we reviewed the derivation of the low-energy effective theory XEFT from its underlying theory QCD for a 
practical description of few-baryon systems at a momentum scale ^ 0{m^). In this section, we summarize results of 
a contact effective theory without pions (EFT(|t), see Sec. VB for references). Its underlying fundamental theory is 
XEFT from which it inherits the QCD-parameter dependence. Because of the stated unresolved problems of XEFT to 
quantify its theoretical uncertainty those errors are also not properly propagated to EFT('^). The set of observables, 
thereby parameterized by couplings of QCE^ and the matching conditions, is smaller than that accessible with 
XEFT which in turn does not cover the measurable parameter space described by QCD. Matching conditions, in 
general, encompass: (i) The point in the QCD parameter space where amplitudes are equated, i.e., the relevant QCD 
interaction parameters {rriq, strong and electromagnetic couplings gs and a); (ii) A set of consistent quantum numbers 
specifying the observables which are expected to be well described in either theory {e.g. , the deuteron, ^a„p, and/or 
the magnetic moment of the triton); (in) A parametrization of the radius of convergence of the effective theory (power 
counting). 

The resultant benefit of a practical theory to express few-nucleon observables in terms of QCD parameters was 
utilizecPi^for the prediction of the three-nucleon spectrum for G (140,200) MeV. In that work, the match to 
EFT(]t) was made to predict three-nucleon scattering observables and excited states of the triton. Prior to Ref. |3T| , 
the ultraviolet RG fixed point in the two-body and the limit cycle in the three-body sector of EFT(^) was conjectured 
at a pion mass close to the physical poiniP^ under the assumption that an amplitude relevant for the description of 
the three-nucleon systerrpj remains constant over the considered range of pion masses which includes the physical 
point from where its value was determined. This assumption was justified a posteriori in Ref. m by solving the 
three-nucleon problem with a XEFT potential and using the obtained pion-mass dependence of to obtain the 
aforementioned first two states of an Efimov spectrum. The limit cycle as a periodic dependence of a coupling 
const anip] is manifest in the neutron-deuteron scattering lengths. In Ref. m the neutron-deuteron scattering length 
was found to diverge in both spin channels at a critical pion mass of ^ 198 MeV with EFT(]t) at LO. This is a 
noteworthy difference to the nucleon-deuteron system at physical m,r, where EFT('/) predicts a limit cycle only for 
the doublet channel and the three-body momentum-dependent counter terms of the quartet channel appear at an 


^ The increasing effect with m-r: of the unknown dependence of the short-range interaction is visualized in Ref. m which in turn will 
also serve as a measure of future rigorous Ci and Di assignments. 

The one-dimensional pion-mass dependence will be augmented once electromagnetism and non-degenerate up and down-quark masses 
are considered in the lattice simulations. 

Originally coined A* in Ref. 00] . 

Only in the doublet channel, this constant can be identified with a momentum-independent six-fermi vertex. In the quartet channel, 
any six-fermi counter term has to be momentum dependent or break (iso)spin symmetry to affect the quartet S-wave amplitude. 
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FIG. 3. (Color online) Pion-mass dependence of the three-nucleon binding energy in the triton channel. For m^< 210 MeV, 
the result^^utilize a match of XEFT to da ta at physical m-n- For m.n-> 300 MeV, LQCD data is shown with box heights 
representing the lattice uncertainty (see Sec. IVB I. The vertical black line marks 771 ^^ ~ 198 MeV at which > oo. 


order expected by dimensional analysis. An ensuing NNLO analysil^ confirmed this behavior of the scattering lengths 
showing traces of a log-periodic behavior in the vicinity of the critical pion mass of about 198 MeV. The dependence 
of the nuclear landscape on ranging from 130 MeV to 806 MeV up to A = 3, the two complementary studies of 
Refs. [31] m] give in combination to the analyses presented below in Sec. VB is shown in Figs. and |^ . 

In the upper row of Fig. the singlet (left) and triplet (right) NN binding energies are shown. The XEFT cal¬ 
culations extrapolate from the physical m,r up to ~ 210 MeV and find a critical pion-mass value where scattering 
lengths in both channels divergq^at m,r ~ 198 MeV. At this value, the deuteron becomes unbound and the singlet-nn 
bound. The LQCD calculations (for references and numerical values see Table [ll|) at 300, 510, and 806 MeV find more 
deeply bound states in both channels. This peculiar situation requires at least one more rriT^ at which diverges 
but allows multiple scenarios for the singlet: (i) tio more critical point and a gradual increase from the shallow state 
(red thick line upper left) to the deeper states; (ii) a intrusion of another bound state at threshold at some critical 
where the former ground state is deeply bound. Both cases indicate the existence of one or more pion masses 
210 MeV < ttItt < 300 MeV where QCD is on the RG trajectory for two-nucleon ultraviolet fixed points and hence a 
three-nucleon limit cycle. To observe signatures of this critical trajectory in LQCD simulations could thus be feasible, 
today. 

What is known about the dependence of the triton binding energy directly (LQCD) and indirectly (yEFT) 
from QCD is shown in Fig. Around the critical point, Ref. [3T] finds the ground state stable (red line) relative 
to the increase observed in LQCD data from 300 to 806 MeV. The above-mentioned tower of excited states which 
forms around is not shown. Instead, we show the manifestations of the unitary fixed point in the three-body 

sector, as calculated in Refs. (SHUT] (XEFT—>'EFT(]t)) and Ref. [33] (LQCD—)'EFT(/)), in the bottom row of Fig. 
For ttItt < 198 MeV, the doublet scattering length ^and (left) is obtained from neutron-deuteron asymptotic states, 
diverges at the critical point to -|-oo. Above the critical pion mass, the deuteron is unbound and is replaced with 


The LECs were deliberately chosen to that end. Singlet and triplet scattering lengths could diverge at different m-^. 
fields could be tuned to realize diverging scattering lengths in all two-nucleon channel was suggested in Ref. 02]. 


That magnetic 
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the then bound nn-singlet state in the definition of The latter rises from —oo. In the neutron-deuteron quartet 

channel (right lower panel), the scattering length goes to +oo at the critical point and cannot be defined in the 
gap between m“‘*' and some below 300 MeV where the state is bound, again. At present, we know of no 
systematic way to connect the XEFT predictions with the LQCD measurements. To understand the emergence and 
disappearance of bound states with LQCD calculations tracing, e.g. , the two-nucleon spectrum over some range 
of m,n., would indicate whether the shallower state at 510 MeV corresponds to the lowered pole at 806 MeV or enters 
at threshold. The statistical uncertainty is the main obstacle to answer the seemingly simple question how reacts 
to an infinitesimal change in ttItt. The slope of as a function of m^n- could be obtained for m^n- > 806 MeV and 
provide valuable information for the extrapolation to smaller pion mass. 

We will see below (discussion around Fig.[^ hints that the nuclear interaction at large is well approximated by a 
SU(4) Wigner-symmetric theory. Assuming that XEFT propagates the explicitly broken SU(2) or SU(3) quark-flavor 
symmetry well into the interaction between nuclei, the XEFT results are not consistent with the LQCD results, which 
use degenerate up and down quark masses. This discrepancy highlights the significance of the symmetry breaking for 
an extrapolation. 

In this section, we summarized results of XEFT and EFT(^) as approximations of QCD. Bound and scattering 
properties of nucleon systems with A < 3 were found sensitive to the QCD-parameter This behavior mediated 
by XEFT was obtained without a systematic power counting and ambiguities due to insufficient data. How a contact 
theory can complement XEFT in the facilitation of scattering calculations while converging in its predictions of the 
three-nucleon spectrum was also shown. Until now, XEFT inherits only its symmetries from QCD. The LECs are 
determined by a match to experimental data at ^ 140 MeV because LQCD cannot be solved for the necessary 
amplitudes for < 300 MeV. However, the variety of nuclear behavior to expect when substituting the scattering 
lengths and effective ranges at physical m^r with LQCD predictions at some m,r ^ 200 MeV where ChPT seems to 
be applicable as matching conditions, was exemplified. 

The analysis of nuclei at larger m^n-, where ChPT is not a reasonable approximation of QCD but where LQCD is 
practical, today, is the subject of the following sections. 


IV. LATTICE QCD FOR MULTI-HADRON SYSTEMS 

To predict nuclear observables rigorously by solving QCD, a chain of systematic approximations has to be employed. 
Each component of this chain is built as an effective field theory and as such provides a prescription to recover the 
underlying theory. The belief in QCD as the relevant theory from which nuclei emerge can be tested by calculating 
observables from correlation functions of the type 

(6) = Z-l J VA^VqVq d{q, q, A^) e" ^ 

The lattice methodology evaluates this path integral and the partition function Z = (I) over fermionic quark fields 
{q,q = and gluon gauge fields (A^) by discretizing Euclidean space time. Thereby, LQCD constitutes the first 

approximation — or chain element in the above terminology — to relate QCD parameters, namely the quark masses 
TOg and the QCD length scale, implicit in the coupling strength to nuclear physics. 

For few hadron systems, in fact, it is the only approximation necessary. For spectral details beyond the ground state, 
larger systems, where O is comprised of > 5 baryons, and/or realistic values of m^, additional EFTs are necessary. 
The various EFTs, ChPT, and EFT(/), for instance, interface through amplitudes with O resembling an operator 
that can be evaluated in both, the underlying and effective theory. 


A. Methodology 

First, the defining parameters of a lattice calculation, roughly taken as the spacing b of the lattice, its space-time 
volume x T, the statistical properties of the gauge-field ensemble, and the to, , gs values affect the accuracy of the 
extraction. Furthermore, an appropriately chosen shape for the operator O increases the accuracy. To access spectra 
of few-nucleon systems, O is composed of hadronic interpolating fields with the generic structure 

^ r^^’"’’°"'’g(ai)... g(a„^) , Oi = {color, flavor, spinor, spatial a;} . (4) 

a 

A hadron is identified with a certain irreducible representation of SU(3) flavor, e.g. , n,p. A, S^’°, 5^ o 8, and 
therefore neither the number of quark fields Uq nor the structure of the tensor F h are uniquely defined by the hadron. 
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Spectroscopic information on hadrons is obtained through the selection of an appropriate operator O, which creates 
a state of defined parity tt, cubic angular momentum I, isospin, strangeness at some initial space-time point and 
annihilates it later. The so called smearing of the quark fields to increase the overlap of the interpolating fields with 
an eigenstate of the QCD Hamiltonian via 


A^{t))q{y, t) 

y 


(5) 


utilizes this freedom (the calculations considered in Table differ in the weight function G and to which parts of O 
the technique is applied). 

All lattice calculations under review, here, employ a source and sink structure for O. States are created with parity 
TT, total angular momentum and Jz, isospin and Iz, strangeness, and baryon number A at the same time to to 
annihilate them at a t > tg- To access the spectrum of a three-baryon system like the triton or Helium-3, for example, 
correlation functions of the form 


CN.N,Nr,riP2,Pi.P3-,t)= E X 

Xl,X2,X3 

X {xiA)N2^ (x2, {X3, t)fVi,/3i (®0, 0)^2,02 (xq, 0)A^3,/33 {xq, 0)) (6) 

are usecP^with spinor indices a, /3. The particle species is understood to be encoded in the Ni’s. The space-time points 
at which the source creates {xqAo = 0) and the sink annihilates {xi,t) hadron(s) with quantum numbers selected 
through r are choices like the projection on a state with defined momenta Pj. Another choice are so called wall sources 
used, e.g. , in Ref. [JS] , which project each quark field of zero momentum out of the vacuum viz. q(t) = q{x,t). 

So far, nuclear observables have been calculated with fields N resembling elements of the baryon octet, only. 
In Ref. [8] , the proton mass is obtained with an A = 1 interpolating field N {x,t) = €abc{u°'’^Cj 5 d^)u^’°'{x,t) 
corresponding to a color (indices abc) singlet which is combined from up (tt) and down (d) quark fields defined at the 
same space-time point but smeared (Eq. ( [^) , and the charge conjugation matrix C acting on the spinor components 
(index a). At the hadronic level in Eq. ™ one chose to define the three sink and source operators at the same 
space-time point. In general, smearing these operators over the spatial coordinates is admissible, too, like shown in 
Eq. § for the individual quarks. 

Back on the hadronic level this dependence on individual quark coordinates is necessary to extract not only spectral 
information, but also wave functions. The standard technique to access both, wave functions '^{Pi 
energy eigenvalues first translates the sink Ns{x,t) = in time before the 

insertion of a complete set of states: 


y-,/r s(ource),si(nk) T \ 

C{{pA '}i=i.....,4i) oc 


(X (A^3i({p/‘}i=i,...,A-l,t)^s({p/}i=l,....A-l,to = 0)) 


E 


r.-Er,t 


2En 


\n){n\ = 0)) 


En ({pz Eji ({pj }i—1; E) ^o(pl,...,x4—11 


— E'nt 


(7) 


For t —>■ oo, the contribution from the lightest hadronic state that couples to source and sink dominates. The overlap 
factors Zn depend, of course, on the sink/source structure and determine the minimal propagation time t for a practical 
extraction of Eq. The residual time dependence in the sink/source overlap with the eigenstates indicates the statistical 
noise which is present through Monte-Carlo sampling in realistic calculations. It is noteworthy, that n = 0 in Eq. 0 
does not need to be the ground state for a given time interval. The state of interest, ground, excited, or continuum 
state, guides the structure chosen for the sink and source. Smearing of quark fields in the nucleon interpolating field, 
for instance, is inspired by the knowledge of the extended nature of the nucleon. For scattering states (see paragraph 
below) in the center-of-mass frame, the correlator with a sink structure which projects hadrons with momenta of equal 
magnitude but opposite direction will be dominated by an excited state over some time interval, because the overlap 
factors Zn in Eq. Q win over the exponential decay. 

As wave functions and energies are calculated in a box of finite size L, scattering phases 6 can be inferred from the 
correlation function. If the energy extracted via Eq. Q corresponds to a scattering state with E = 2k'^/m^, Liischer’s 
formula (original work in Refs. [3S1137] , for boosted systems see Ref. [JB] ) 


A:coti5 




- 47rAj 


( 8 ) 
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TABLE I. Parameters of the LQCD measurements reviewed in this work. Few-nucleon systems are calculated on a L® x T 
lattice with spacing b using three or two degenerate quark masses with resultant pion (m^) and nucleon {rritq) masses in an 
ensemble comprising A^^fg gauge configurations. 


QCD version 

L [fm] 

T [fm] 

b [fm] 

m.,r [MeV] 

ruN [MeV] 

Ncfg 

SU(3)PIS3 

3.4 ^ 6.7 

6.7^ 9.0 

0.15 

806 

1634 

> 1905 

SU(2)P 

2.9 

5.8 

0.09 

701 

1583 

390 

SU(2)PI 

2.9 ^ 5.8 

4.3 ^ 5.8 

0.09 

510 

1320 

200 

SU(2)P 

2.5 

4.0 

0.13 

354,493,593 

- 

490 660 

SU(2)P 

2.8 ^ 5.6 

7.5 ^ 11.2 

0.12 

450 

1226 

> 1000 

SU(2)P 

4.3 & 5.8 

4.3 & 5.8 

0.09 

300 

1050 

400 & 160 

gU(3jp[II] 

3.9 

3.9 

0.12 

469 ^ 1171 

1161 ^ 2274 

720 ^ 420 


with a sum over all integer three-vectors j with magnitude smaller than the lattice-momentum cutoff A, relates the 
energies to phase shifts. 

The algorithm employed by the NPLQCD collaboration to calculate multi-baryon correlation functions summarizes 
this section. 

First (hadron level), an A-baryon operator is defined via Eq. Q with flavors taken as elements of the baryon octet. 
All baryons comprising the source in Refs. [Hill] are defined at a single lattice site, i.e., the sum in Eq. Q excludes 
spatial degrees of freedom, and the fields are created at the same Xq. The sink is chosen to project out a state 
with defined momentum for each individual baryon, i.e., the sum in Eq. runs over all lattice sites and the factor 
contains a plane wave for each baryon. In principle, the ensuing momentum dependence of the overlap 
function Eq. Q provides a rigorous assessment of the typical baryon momenta within a nucleus. Measurements like 
this are desirable for the construction of effective theories (see Sec. VB) as they could provide their typical scales. 

Second (quark level), the baryon operators are substituted with smearecP 13 (smeared and poiniP) interpolating 
operators, again defined at a single space-time point. 

Erom correlation functions constructed in this way, properties of nuclei with A < 4 have been calculated at the 
SU(3) flavor-symmetry point and for degenerate up and down quarks with the strange at its physical mass (SU(2)). 
We compare the parameters of those calculations in Table |T] using the lattice spacing b, the temporal and spatial 
lattice size T x L^, the number of gauge configurations iVcfg, and the QCD version, i.e., SU(3) or SU(2) flavor, as 
standard. The selected calculations supersede the pioneering studie^^^Ml-vyhich rely on the uncontrolled quenched 
approximation. In contrast, the analyses under review in Sec. |IVB] use controlled approximations, only, in the sense 
that the usage of more appropriate parameters as the ones given in Table |l] must yield results consistent with the old 
predictions within error bars. 


B. Data 

The data selected for this review represents the most advanced LQCD extractions (judged by the standard pa¬ 
rameter set listed in Table |l]) of spectra in various two, three, and four-nucleon channels which employ controlled 
approximations, only. The two extractions of Refs. [ 81 E 2 ] use the same mass for all three light quarks, set such that 
^ 806 MeV. The exploratory calculation three-baryon systems in Ref. [H] , Ref. [7] (510 MeV), Ref. [5] (300 MeV), 
Ref. |3D] (350 —590 MeV), and Ref. [5] (450 MeV) use the physical value for the strange-quark mass but degenerate 
up and down-quark masses corresponding to the respective At the SU(3) symmetric point, data over a wide 
range of states, e.g. , the H-dibaryon (AA), the hyper triton (^H), or hyper Helium-4 (AHe), was extracted. No 
comprehensive theory has been devised for few-baryon systems comprised of all elements of the A = 1 octet based on 
this data and we review the strangeness s = 0 sector, only. We begin with the two-body sector where we compare 
features of the nuclear bound and scattering systems at physical with the LQCD data. The three and four-nucleon 
sectors are discussed in parallel, given their correlation. 

Two nucleons — Eventually, LQCD must produce the two small (relative to Aqcd 47r/^ ~ ~ 1 GeV set by 

5s(Aqcd) 1) scales characteristic for nuclear physics in their observed unnatural ratio, namely, a momentum scale as¬ 
sociated with the poles of the low-energy two-nucleon scattering amplitude, (^a„p)“^ ~ = 73 ^ 45 MeV 

in the S-wave spin-triplet and (^Unp)”^ ~ 7 ]^ ~ —8 MeV in the S-wave spin-singlet channel, and another related to 
the typical range of the nuclear interaction ^ Tn~^. In addition to this unnaturally large value, m 7 r /73 ~ 3, only the 
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TABLE II. Few-nucleon data from physical and computational experiments. LQCD uncertainties are upper bounds of the 
quoted values. Binding energies, are given in MeV, 73 is the binding momentum of the two-nucleon triplet state, 

and B 2 are the binding energies corresponding to the poles of the effective-range amplitude (see Eq. parameterized by the 
central values of the scattering length ^a„p and the effective range r^. Data at niTr ~ 806 MeV has been updated^. 



TABLE III. Two-nucleon scattering lengths, (in fm), and effective ranges, r 2 s-i-i (in fm), above the inverse physical 

and unphysical pion masses (in fm). The data sources are listed in Table [I] 



140 

353 

450 

493 

593 

806 

(Inp 

-23.75 

0.63 ±0.50 

20 ±60 

0.65 ±0.18 

0.0 ±0.5 

2.3 ±0.46 

ri 

2.75 

- 

3.0 ± 1.0 

- 

- 

1.1 ±0.14 

CO 

73 

5.42 

0.63 ±0.74 

-11 ±55 

0.41 ±0.28 

-0.2 ± 1.3 

1.8 ±0.31 

73 

1.74 

- 

3.4 ± 1.8 

- 

- 

0.91 ±0.14 


1.4 

0.56 

0.44 

0.40 

0.33 

0.25 


triton binding momentum has to be known for a predictive theory whose breakdown with the nucleon number A is 
still unknown. 

Associating scales analogously at heavier this parameter, i.e., in general the scale of the system compared to the 
range of the interaction, and here in particular 771 , 7 / 7 , is still large (Table [HI 6*^ row) and thus suggests an unnatural 
EFT. At physical m,,, the scattering length and effective range provide another pair of scales that approximate the 
low-energy two-nucleon spectrum equally well as a consequence of the pion dominating (ex- or implicitly) the nuclear 
interaction up to ~ 100 MeV, although it cannot account for the emergence of the unnatural size of the system. Both 
measures, and a/r, are therefore of similar size. This is also found at tti,, 450 MeV, suggesting a similar 

analytic structure of the amplitude. The first fully-dynamical calculation of NN scattering parameters in Ref. [ 30 ] 
extracts scattering lengths, only. These scattering lengths are ^ (see Table III) and thus indicate natural NN 
systems at 350, 490, and 590 MeV pion masses in both spin channels. Here it is assumed that the interaction range is 
about the inverse pion mass. This assumption fails as shown in Ref. |52j , where ~ 806 MeV does not approximate 
the effective interaction range well. By taking the lattice data for a and r as scales, a more natural theory in which 
these parameters are of the same size is implied. Comparing the respective ratio with scales given by the system’s 
binding momentum and m,,, yields a larger ratio. 

A graphical summary of available scattering-length measurements as compiled in Table jlflj is shown in Fig.[^ There, 


we display the experimental 
larger pion mass 


np 


^a„p: filled red. 


(filled black circle) and ^Unp (empty black circle) together with the lattice data at 
Gnp'. empty blue). The continuation (transparent areas for tti,, > 350 MeV) of 
the allowed region of scattering-lengths (opaque blue for triplet and opaque red for singlet) as extrapolated^ using 
the BBSvK power counting of Ref. m , is a speculation which is motivated by the log-periodic running of the 
three-nucleon momentum-independent interaction with a regulator parameter (see next paragraph and Ref. |57j ) . As 
shown in Fig.j^ the NPLQCD data indicates a critical pion mass of about 440 MeV in addition to the one founcj^at 
198 MeV at which NN scattering lengths diverge. The limit-cycle assumption of a periodic behavior of two-nucleon 
scattering lengths, with the pion mass instead of a cutoff parameter, leads to a nuclear interaction at a larger pion 
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FIG. 4. (Color online) Pion-mass dependence of the singlet (filled circles) and triplet (empty circle) neutron-proton scattering 
lengths a. Allowed regions for a for < 350 MeV (opaqne) were derived^ with KSW power conntin^^. The transparent 
regions are speculations inspired by the limit cycle as observed in the running of the three-body contact interaction. 


mass identical to the one found in nature. Real nuclej^ could thus be calculated from QCD on the lattice at large 
pion masses, avoiding the uncertainty associated with light quark masses. From Fig. one would naively expect this 
copy of the real nuclear two-nucleon world at about 850 MeV. 

We turn the discussion to data sets which include effective ranges. Scales are related to the analytic structure of an 
amplitude, i.e., the position of a pole and the radius of convergence of some expansion around that pole. In case of the 
effective range expansion at physical m^n-, 7 ^ a~^ and r ~ respectively. In that sense, 7 ~ ^ k±{a,r) 

(see poles of Eq. ^ below) defines a scale even at ~ 806 MeV. But 7 ~ ^ k± (a, r = 0) and TO7r/7 ^ a/r 

(Table [ll| imply additional non-analyticities within an 771 ^^ radius around 7 beside the existing bound-state pole. 

To parameterize the scattering amplitude using effective-range theory (ERT) as 


Tik) 


47r 1 

Ton k cot 6 — ik 


4tt 1 

-b ^k^ + 0{k^) - ik 


(9) 


is justified through a pole consistent with the non-perturbative binding energy (compare the 1''*(2”‘*) and 7“'(8*‘*) row 

possibly two bound 
is beyond the range of 


in Table. . The amplitude Eq. 
states. At physical m,r, k'^/mr. 


has two poles at momenta k± = y(l ± y/l — ‘irja), 
2.2 MeV, the binding energy of the deuteron, while 


validity of the ERT and therefore not in disagreement with the missing experimental evidence for such a deep state. 
At TOtt 806 MeV, 2r/a ~ 1 in both S-wave channels (see also discussion in Sec. VBI which implies relatively closely 


14 


Of course, the role of the nucleon mass, and the effective range has to be taken into account. 
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FIG. 5. (Color online) Dependence of two-nucleon binding energies as poles of the ERT scattering amplitude parameterizecP^at 
~ 806 MeV via ^a„p = 2.33 ± 0.46 fm, r(^5'o) = 1.13 ± 0.14 fm, and = 1.83 ± 0.31 fm, r(^Si) = 0.91 ± 0.14 fm. r — a 
pairs from the gray lattice area in the xy plane allow virtual or resonant states, r — a pairs defining an e nerg y (purple 
patch of the surface) consistent with the value extracted with a “bound-state” source (see in Tabled and a gap to 

^2 ~ 30 MeV (light gray patch) are marked in green. Binding energies with Sj — < 30 MeV are shown as red surfaces. 

The two surfaces touch at a line which marks the double pole of the amplitude. 


spaced bound states. Considering the uncertainty in the lattice extraction of the scattering length and effective range 
yields via Eq. @ the following scenarios of the two-nucleon spectrum: (i) two shallow states for r/a = (2-be) ^ (dark 
and light red areas on the energy surface in Fig.[^, (ii) one shallow state and a spurious deep state if r/a —>■ 0.3, z.e., 
the analog to physical (gray energy surface and its green projection onto the a — r plane), and (in) no bound 
states but resonances for rja ^ 0.6 (hatched gray area in a — r plane). The limits correspond to the smallest and 
largest ratios consistent with the uncertainties. 

Scenarios (i) and (ii) are also consistent with the spectral measurement of = 22.3±5 MeV albeit only one bound 
state has been isolated. To discriminate a second, almost degenerate state is at present impossible as the individual 
extraction of energy levels from the effective-mass plots suffers from significant statistical noise. A ratio measurement 
between a singlet and a triplet two-nucleon bound state correlator will be less affected by that noise and could hint 
towards a second bound state. In Sec. VB we comment on the implication of this peculiar analytic structure for the 
construction of an EFT. 

Three and four nucleons — For heavier systems, LQCD provides, at present, only data on bound states, until 
numerical precision will suffice to identify signals corresponding to two-fragment continuum states thus enabling the 
calculation of, e.g. , nucleon-deuteron phase shifts via a generalizatioiJ^of Eq. 

At physical the existence of bound states in the three-nucleon system is a universal consequenc^^of the 
proximity of the two-nucleon system to the unitary limit —>■ oo. The fact that there is only one such state in 


the channel and its exact binding energy are then peculiar to the three-nucleon system and not fixed by the 

two-body scattering length and effective range. Furthermore, there is no physical bound state in the | channel. In 
contrast to this independence of of the two-body a’s and r’s, the three and four-nucleon spectra are intertwined, 
and the emergence of two four-body states with every three-body state, as known from unitary bosonic system^^, 
seems to generalize to the nuclear problenJsni . Whether or not the Coulomb repulsion transforms one of the two 
bound J'^ = O"*" states to a shallow resonance and the other to the a-particle will be a test for EFT(^) and its 


















Matching effective few-nucleon theories to QCD 


14 


TABLE IV. Binding energies (in MeV) of the shallowest two, three, and four-nucleon state relative to its lowest break-up 
threshold, i.e., A21 = B 2 CS 0 ) {B 2 {^Si) at ~ 140 MeV), A32 = B3(^Si/2) - B 2 {^Si), and A43 = B^CSo) - BsCSi/i). 
Uncertainties (insignificant at rtir: ~ 140 MeV) were combined in quadrature. References to LQCD data at unphysical are 
the same as in Table H) 


rUn 

A 21 

A 32 

A 43 

A 43 /A 32 

140 

2.2 

6.3 

20 

3.2 

300 

8.5 ±2.3 

7.2 ± 17 

25 ±41 

3.5 ± 14 

510 

7.4 ± 1.4 

8.8 ±5.8 

23± 19 

2.6 ±3.8 

806 

16 ± 5.4 

33± 11 

44± 24 

1.4 ± 1.2 


treatment of electromagnetism in few-body systems. 

The LQCD calculations at higher m,r identify single bound states in the triton and a channels, too. There is no 
data on other three or four-nucleon channels. In the strangeness s = —1 sector, the NPLQCD collaboratioiP finds a 
bound hyper triton. At the considered SU(3) symmetric point, this state hints at a corresponding bound s = 0 
quartet state. For s = —2, a second four-body 0"*" bound state is found, reminiscent of the aforementioned universal 
tetramer pair associated with each three-body state in the two-body unitary limit. Given the degenerate and 
i'*' three-body, and a and aaKo four-body states, and the above conjectured two-level structure of the two-nucleon 

system, LQCD investigations of possible three and additional 0+ four-nucleon bound states could reveal striking 
differences of few-nucleon systems at larger pion masses. 

Beside the sheer existence, the relation of binding energies, i.e., thresholds, and binding energies per nucleon in 
nuclei of different A are relevant for an effective description of heavier systems. At physical m.^-, the triton is bound 
by ~ 6.5 MeV relative to the deuteron-neutron threshold and interpreted as shallow, corresponding to an Efimov 
state. The a—particle is bound by ~ 20 MeV relative to triton-proton and considered a universal feature. Some 
hints to whether or not this interplay between two, three, and four-nucleon states persists at higher m,^. are found in 
the relative threshold positions, compiled in Table |IV[ For all ttItt > 140 MeV, the triton is closer to the deuteron- 
nucleon threshold than the a—particle is to the noninteracting triton-proton system. The three-nucleon (A 32 ) and 
two-nucleon break-up energies are of the same order and decrease from 806 to 510 MeV. At ttItt ~ 300 MeV, a bound 
or an unbound three-nucleon system are within uncertainty limits. 

The gap between the measured bound state energy of a and its lowest break-up threshold A43 is larger than A32 
but of the same size or smaller than the scale set by i?^. The ratio A43/A32 decreasing with ttItt and A43 increasing 
simultaneously, is reminiscent of the scattering-length dependence of the two and three-body systems for a —>■ 00 . 
With decreasing m,r, the triton approaches threshold. The discussion below on the neutron-deuteron scattering length 
suggests a diverging three-nucleon amplitude at zero energy and thereby the analog of a limit cycle in the four-nucleon, 
three-body deuteron-nucleon-nucleon system. 


We summarize this section at the beginning of Sec. VB which reviews a theory trying to describe the data consis¬ 
tently. 


V. NUCLEAR THEORIES 


The interpretation of the above results as a data base for a theoretical analysis is analogous to the way early 
experiments on nuclei initiated theoretical nuclear physics. A theoretical analysis of lattice measurements is justified 
if computational resources are not expected to be available in the near future for the nuclear properties of interest. 
Two approaches to a systematic understanding of nuclear lattice data are available. In effect, they generalize the 
concept presented in Sec. |III| by matching a nuclear contact theory to LQCD amplitudes and the application of the 
ensuing theory to few-nucleon systems. The two methods differ in the matching condition and the contact theory. 
How exactly those differences lead to inconsistent postdictions is not knowrp^ Thus, we deem a brief summary of 
their respective technique and basic assumptions as useful. 


For recent work on the assessment of the sensitivity of few-hadron LQCD results on the source structure see Ref. m ■ 
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A. Matching wave functions 


Analogous to an approach taken for the description of Kaon decaj®, a connection between QCD four-point 
correlation functions and non-relativistic nuclear potentials was made (the original work is Ref. [45j , for a review see 
Ref. |M1 )• While matching as reviewed between AEFT and experiment or EFT(]t) identified amplitudes, here, two 
theories are matched through a set of wave functions. The underlying theory QCD defines a wave function which can 
be extracted with the lattice technology. This function satisfies a non-relativistic Schrodinger equation if the potential 
is chosen appropriately. The unknown is thus the potential while the wave function is input. 

More specifically, a relativistic Nambu-Bethe-Salpeter wave function (introduced in Ref. [5S] for single hadrons) 
is extracted with a correlation function as in Eq. ([^. Ref. [IS] , in particular, uses a so-called wall source for 
A^source({p/°“‘^'^®}i=i,...,A-i, to = 0), 1 . 6 ., the local quark fields in the interpolating operators are projected onto zero 
momentum. The source specifically creates a two-nucleon state with defined parity, isospin I and total angular 
momentum J. At the sink, the two-nucleons are annihilated with the same quantum numbers but at a selected 
distance x apart: 


C{x, t-to) = {Nijj{r, t)N 2 j,j{r + x, t)lVi,/,j(to)A^ 2 ,/,j(to)) 


( 10 ) 


We iterate the aforementioned ambiguities in the calculation of this correlation function: the structure of the inter¬ 
polating fields at the quark and hadron level. 

Two methods have been used to extract a two-nucleon potential from the QCD Greens function C{x,t — to) in 
Eq. (10). The standard used to predict masses and ground-state energies of composite objects in LQCD uses the 
spectral decomposition of the correlation function as given above in Eq. ([^ which reads in the case of interest here. 


C{x, t-to) =Y,Zu- ^'(a:; fe„, I, J) ■ 

n 

Z{ko, I, J) ■ 4'(a;; ko, I, J) ■ e" ^ 


( 11 ) 


For large enough t — to the ground state can be isolated if the statistical noise is sufficiently reduced. To reduce this 
noise is a challenge for for the identification of energies from effective-mass plots, in general, and can be achieved 
through an increased m.^. It is problematic for the extraction of wave functions, in particular, as shown with a toy 
model in Ref. [66j , which exemplifies how admixtures to the wave function 4'(a;; fco) J) stemming presumably from 
other states in the spectrum modify a potential derived from it. That model, specifically, found the hard core, relative 
to the outer tail, of the potential to be highly sensitive to the coupling to other channels. Thus the amount of ground 
state saturation in Eq. 0 controls the range of applicability of the potential — the stronger the admixture, the 
smaller the energy at which the potential becomes useless. 

An alternative extraction method which overcomes the saturation issue was devised in Ref. (Ml (see Ref. [57] for a 


comparison). It cleverly exploits the fact that each wave function in the spectral sum Eq. (111 satisfies a Schrodinger 
equation. Therefore, one can use a rescaled correlation function R{x,t) = (^(a:, t)/(e“™”')^ to define a potential U 
which, in general, acts between two-nucleon channels and thus shall be, like C, understood as a matrix, via 


d 


-Ho-^ 


1 


dt drriM dt"^ 


Ra{x,t) = J d^yUaa'{x,y)Ra'{y,t) 


( 12 ) 


The same potential can be constructed from the ground-state wave function obtained via Eq. 0 from the stationary 
Schrodinger equation 


{-ffo + ^} ^'a(a:) = ^ J d^yUaa'{x,y)'i>a'{y) ■ 


(13) 


A non-local, velocity-dependent ansatz 

U{x,y) = V{x, V)5^‘^\x -y) = Vo + Kcri • era -b VtSi 2 + • S -I- C>(V^) 

was choseii^for the nuclear potential in a given channel. The expansion shown on the right-hand side defines 
I a; I-dependent coefficient functions which are constructed order by order inverting either Eq. (12) or Eq. (13). 


(14) 


Nuclear observables in the strangeness s = 0 sector predicted with potentials as defined above relying on unquenched 
lattice QCD correlation functions (see Tableinclude the scattering length in the two-nucleon channel with c.m. 
momentum k derived in the limit 


(tott = 701 MeV) = lim k ^tani5(A:) = 1.6 ± 1.1 fm , 


fe-j-O 


(15) 
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FIG. 6. (Color online) Pion-mass dependence of momentum scales relevant for nuclear low-energy physics. The green-shaded 
boxes represent the range of binding momenta y/2m^^~BjA for nuclei with A = 2,3,4 and binding energy B as given in Table [n[ 
The pion-mass (black) and A (gray, see footnote 161 scales become equal for m.n- ~ 806 MeV. 


which was obtained from the central part in Eq. (14|. Tensor and spin-orbit potentials have also been derived^ at 
up to 1.1 GeV and applied in variational calculation!^. Neither two nor three-nucleon bound states but indications 
for a shallow four-nucleon state, whose binding energy increases from ~ 0.8 MeV to ^ 5.1 MeV with m,r decreasing 
over the considered range, have been found. Employing the same potentials via the Brueckner-Hartree-Fock method 
to and ^°Ca nuclei yields those nuclei bound consistent with the 5.1 MeV for the a-particle. 

Euture work deriving a three-nucleon potential, as pioneered at "^1.1 GeV in Ref. |69| , will tell whether or 
not the above expansion of the nuclear interaction receives a significant contribution from such a force relative to 
higher-orders of the velocity expansion of the two-nucleon potential. In the following section VB[ we will present an 
analysis of the role of three-nucleon force and its dependence in an EET framework. 


B. Matching correlation functions 

A theory for few-nucleon systems which explains the LQCD data at niTr ~ 510 MeV and 806 MeV (for lattice 

data see Refs. [ZlIH] , and for the EFT Refs. ^03] ) resembles the effective-field-theory approach to low-energy 
few-nucleon systems in the physical worlcP^ Compared to the method introduced in the previous section, it 

matches to QCD via observable binding energies instead of wave functions. The matching conditions, as part of the 
definition of an EFT, replace real-world data with LQCD predictions for the input — we are exploring the domain of 
the transparent sheets in Fig. 

Nucleons as eigenstates of QCD are canonically defined as an isospin doublet belonging to the lowest-mass baryon 
SU(3) octet. The first amendment to the most general SU(3) invariant Lagrangian, constructed solely from this 
octet, couples it, first, to the lowest-mass meson octet and second, to the lowest-mass baryon decupl et. A comparison 
between (i) the nucleon mass Ton, a scale associated with the excitation of a nucleon to a ^2m^(m\ — ttin), 
(Hi) the mass of the pion m^^., and (iv) the binding momenta of an A-body nucleus, ^J2m^B jA gives an indication 
whether an approximation of those couplings, specifically, pion-nucleon, A-nucleon, and A-pion-nucleon, with contact 

















Matching effective few-nucleon theories to QCD 


17 


0.8 



400 800 1000 1200 1600 

A [MeV] 

FIG. 7. (Color online) Dimensionless couplings of the lowest-order two-body operators of EFT(^) at m,r ~ 510 MeV (hatched), 
niTr ~ 806 MeV (solid), and physical m,r (red solid line). The upper (lower) bands show the regulator dependence of the SU(4) 
symmetric (asymmetric) LEG defined in Eq. ( |l6| . The band widths represent uncertainty in the data used for the 

calibration. 


interactions amongst nucleons might be useful over some energy range for nuclear amplitudes. 

These scales are compiled in Fig. |^for the different pion masses and led to the ansatJ^of a nuclear contact theory 
analog to the established EFT(^) at physical m,r. In essence, this analogy utilizes: (i) the much smaller typical binding 
momenta of nuclei relative to the nucleon mass to justify a non-relativistic treatment implying small-momentum 
Lorentz symmetry; (ii) with the typical momenta in bound systems — identified with Q — up to A = 4 being smaller 
than the lightest meson and lowest baryon excitatior|^ TOn — WA — collectively denoted M — an (iso)spin 1/2 
nucleon field {N = {n,p)) suffices as sole degree of freedom; (Hi) the renormalized two-body LECs are assumed to be 
of order C 2 n = hence only iterations of the two u = 0 zero-derivative interactions are of the same order. 

The infinite iteration introduces the two-nucleon bound state poles; (iv) the three-nucleon momentum-independent 
interaction in the ^S' 1/2 (triton) channel is considered at the same order as the momentum-independent two-nucleon 
interactions in the and ^Si channels. This unnatural enhancement and the non-perturbative treatment assumes 
that the leading three-body interaction follows a limit cycle analogous to physical in Ref. m . This approach was 
demanded by the cutoff dependence of the triton ground state with only the leading-order two-nucleon interaction. 
Two observations support the underlying assumption of a limit cycle: first, the parabolic increase of the triton ground 
state energy as a function of the cutoff in the absence of a three-body force, and second, the appearance of an 
additional state at some critical cutoff value. 

The Lagrange densitj|^ defining the theory together with the size estimates for the LECs in an operator basis that 
splits SU(4) symmetric (Cq ) and asymmetric (Cg) components reads 

£ = N^ (zdt + Vy(2 m«))N-(N^Nf+ C;[ (N^a-Nf^ +Di(N^Nf . (16) 


The scale in Fig. in contrast, considers the effect of intermediate states under natural assumptions for the LECs. Under these 
assumptions, m-^ remains the lowest threshold setting the convergence rate of the theory. 

Space-time and field coordinates are rescaled such that dt and V^/(27nN) are of the same order as in NRQCD (see comment below). 
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TABLE V. Leading-order EFT(/) post (m,r ~ 140 MeV, B^) and predictions^ for the quartet and doublet neutron-deuteron 
scattering lengths *‘a„d and and for the a-particle binding energy at three pion masses. 


mn [MeV] 

140 

510 

806 

^Und [fm] 

5.5 ± 1.3 

2.3 ± 1.3 

1.6 ± 1.3 

^a„d [fm] 

0.61 ±0.50 

2.2 ±2.1 

0.62 ± 1.0 

R„ [MeV] 

24.9 ±4.3 

35 ±22 

94 ±45 


The EFT and power-counting scheme implied by (i-iv) was implementecSI^ as a cutojf EFT in coordinate space. 
The employed renormalization comprised (i) the regularization of contact interactions with Gaussian functions on the 

relative coordinate S{r) —)■ g^3/2 e~^^^, and (ii) the calibration of the A-dependent LECs to B^, and B^. The 
necessary calculations numerically solved the appropriate two and three-body Schrodinger equations. 

An approximate Wigner SU(4) symmetry of nuclear interactions at m,r 806 MeV has been noted in light of 
the degenerate (within lattice uncertainty) two-nucleon and ground-state energies in Ref. [8] and scattering 
length to effective range ratios a/r ^ 2 in Ref. in the respective channels. This independence of the nuclear force 
w.r.t. SU(4) spin-isospin transformation^^ SN = translated into the LECs as shown in Fig. While the 

relatively small SU(4) asymmetric LEG C^can be understood from the large scattering lengths in both channel^^at 
physical (lower red solid line in Fig. the degeneracy ^a„p ~ ~ 2 fm is explicit at ttItt ~ 806 MeV. At 

physical SU(4) invariance results from being close to the unitary limit, while at larger m,n., it seems to emerge as 
a unique feature of QCD! 

The power counting is justified via RG invariance. Namely, for a cutoff EFT, predictions must converge at every 
order if the RG-flow parameter A —>■ oo. Furthermore, cutoff dependence at a given order can be eliminated at some 
higher, not necessarily the next, order, where an LEG with the right scaling counters the dominating A dependence 
which is not eliminated by the lower-order LECj^ At present, leading-order predictions with A = 2 — 8 fm“^ are 
available with no sign of inconsistencies in the power counting. Specifically, the a binding energy the nucleon- 
deuteron scattering lengths in doublet and quartet were analyzed in that light. The former as a bound-state 
observable that can be benchmarked with the lattice data, and the latter as scattering properties which constitute 
predictions, LQCD wilPmeasure given the resources to increase its numerical accuracy. The predicted observables 
are compiled in Table [V] and also included as gray columns in Fig. for B^. In contrast to physical m.^., where 
the uncertainty of real-world experiments is insignificant relative to the absolute observed value, EFT('^) applied to 
lattice data has to propagate the uncertainty in the amplitudes used to renormalize the theory, namely B^^^ B^, and 


Rj. To justify the power counting, A-variation suffices but was showiPl^to converge slowly. Calculations within 
the dibaryon formalism, which is more flexible in its regulator, could probe the sensitivity of nuclear systems at high 
quark masses to short-distance structure more comprehensively. 

A different power counting altogether is required for the case of degenerate {r/a = 1/2) or closely-spaced 
{r/a = 1/2 ± e) two-nucleon bound states in the same spin channel. In the case of r/a = 1/2, the effective-range 
amplitude Eq. (|^ has a double pole which cannot result from an iteration of momentum-independent terms as in 
Eq. (161. The breakdown of a contact as an approximation of a Yukawa theorjj^with the latter sustaining a shallow 
and an excited state was demonstrated in Ref. m . The similarities indicated in Ref. m of the presumed two-nucleon 
LQCD spectrum with more than one shallow state to non-relativistic QCD (NRQCD) become more detailed in the 
description of heavy quarkonium in NRQCEP^, i.e., bound states of heavy (mass M) quarks and antiquarks with 
typical binding momenta 1 / MB and splitting between radial excitations of order 2R (note the splitting between R^ 
and B 2 in Table B. The EFT developed in Ref. I7H] from full QCD for quarkonium also considers dt and V^/(2M) 
of the same order( compare LO EFT(-/-) in Eq. (16) which uses the same counting for a shallow nucleon state), is 
approximately quark-spin independent, but retains the coupling to Coulomb and transverse gluon fields. The latter 
as a systematic way to incorporate the model scalar field as used in Ref. m from QCD could serve as an alternative 
to the EFT(-/) approach if the second state is found on the lattice. Another way to implement an excited state in 
an EFT considers it as an additional degree of freedom. This leads to the so called dibaryon formulation of EFT(-/), 
which introduces the and the ^Si bound state as degrees of freedom appropriately coupled to the single-nucleon 
field (application in nuclear physics: Refs. [73 jTS] , general quasi-particle mechanism: Ref. izg ). 

Returning to the discussion of the results obtained with the EFT appropriate for the data as available at present, 


G {1,2,3, 4}, a. parameterizes an infinitesimal SU(4) transformation with SU(2) generators acting on (iso)spin 

degrees of freedom. 

We understand A as a placeholder for any cutoff introduced to regulate the theory. It is common, e.g. , in solving the three-body 
problem with the dibaryon formalism to use a combination of cutoff and dimensional regularization. In that case, the limit A ^ oo 
must be taken in both regularization schemes. 

This is to be understood as any mechanism that produces a Yukawa potential, e.g. , the coupling of the nucleon to a scalar field. 
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For each system, the physical (red) value is compared with LQCD measurements at ~ 300 (blue), 450 (hatched purple), 510 
(green), and 806 MeV (purple). Postdictions for of EFT(^) are in gray. The black box represents XEFT prediction^^Uajj;;! 
its height the considered niTr range of (190,210) MeV. The other box heights resemble the lattice uncertainties. 


i.e., with single bound two-nucleon states in both channels, the analysis in Ref. [43] concludes with the calculation of 
Phillips and Tjon correlations. The former is between the triton binding energy and the ^5'i/2 neutron-deuteron 
scattering length, the latter between and B^. While the one-parameter correlations exist by constructiorj^ the 
relatively small sensitivity w.r.t. cutoff variation suggests that the above mentioned entanglement of three and four- 
nucleon spectra persists. The Phillips line, in addition, was calculated down to values for B^ close to the deuteron 
breakup threshold. The observed divergence demonstrated the resemblance of the three and two-nucleon systems, 
where in the latter, the diverging scattering length for the bound state energy approaching zero is well understood. 

Addressing the breakdown of the conceived EFT with the particle number, the ®He and ®Li ground states were 
analyzed^. ®He was found bound for small and unbound for larger cutoffs implying an improper renormalization of 
which signatures were also found at physical m,r in Ref. [80] . For ®Li, results are limited to a single cutoff value and 
bind the system with the same energy per nucleon as in a. Both, the five and six-nucleon system as characteristic 
features of nuclei pose challenges to any nuclear theory: the former for its inability to sustain a bound isotope, the 
latter for its halo and borromean nature. To relate those properties to either two-body ERE and triton parameters or a 
larger set which includes information about the short-distance structure, e.g. , P-wave NN phase shifts or four-nucleon 
resonance parameters, poses a problem which is unresolved at physical m,^., too. The analysis of the dependence and 
evolution of 5 and 6-body peculiarities on drives the development of EFTs for larger, physical nuclei at those 
crucial nucleon numbers. 
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TABLE VI. LQCD predictions for nuclear magnetic moment^^/r, polarizabilitie, scattering length a, and effective 
rangd^r at rUn ~ 806 MeV. Magnetic moments are compared with physical value^^. Uncertainties are naive sums of the 
combined energy shift extraction and magnetic-field-dependence fit uncertainties plus finite-volume effects. In the gray cells, 
results follow naive shell-model expectations. 


State 

a [fm] 

r [fm] 

P ‘806 [nNM] 

pi4o [nNM] 

jg(MO) 

n 

- 

- 

-1.981 ±0.023 

-1.91 

1.253 ±0.12 

p 

- 

- 

±3.119 ±0.097 

±2.80 

5.22 ±0.89 

nn 

2.33 ±0.46 

1.13± 0.14 

- 

- 

1.872 ±0.20 

pp 

2.33 ±0.46 

1.13± 0.14 

- 

- 

5.31 ±2.8 

d 

1.83 ±0.31 

0.91 ± 0.14 

±1.218 ± 0.125 

±0.857 

4.4 ± 1.8 


- 

- 

±3.56 ±0.23 

±2.98 

2.6 ± 1.8 

®He 

- 

- 

-2.29 ±0.15 

-2.13 

5.4 ± 2.4 

a 

- 

- 

- 

- 

3.4 ± 2.2 


VI. ELECTROMAGNETISM 


While for single hadrons the effect of the U(I) gauge symmetry was calculated with the lattice methodology {e.g. , 
the proton-neutron mass splittin^^), the application to systems with A > 1 remains impractical. What became 
feasible is the analysis of the response of nuclei with A = 1,2, 3,4 to a static external electromagnetic fielcP^in the 
form of magnetic moment^^and polarizabilities, and the capture process of a neutron by a proton via the emission 
of a photon np —^ d'j^. In this section, we review the interplay between LQCD data and contact theories in the 
extraction of those observables. 

Data — Here, as sketched above in Sec. |IV A[ correlation functions of operators with non-zero overlap with the 
states of interest provide the spectral information. In Ref. |82j , a uniform magnetic background is considered with 
the approximation of zero-sea-quark electric charges. The approximation allows to recycle the functional determinant 
which results from the analytical integration over the quark degrees of freedom in Eq. ([^. It includes the U(l) gauge 
field in addition to the gluons only in the propagators for the valence quarks which stem from the contractions of the 
fields in the source and sink operators O. In comparison to the spectral and scattering calculations (first row Table |^, 
a smaller lattice, x T ^ (3.5 fm)^ x 5.3 fm, but the same m,^. ~ 806 MeV and three mass-degenerate quark flavors 
(SU(3)) with a value corresponding to the physical strange-quark mass were used. 

The energy levels were obtained in this partially-quenched approximation from effective-mass plots as a function of 
the magnetic field strength |B|. By relating these levels to the expected energy eigenvalues for a 
hadron h G {neutron, proton, deuteron, nn,pp, ^H, with charge Qh, a spin j < I, J ~ e^, and zero mo¬ 

mentum occupying the Landau level 

Eh-jA^) = Ml + i2nL + l)\QheB\ 

-M, • B - 27r^r°)|Bp - + 0(|B|4) , (17) 

the magnetic moment magnetic scalar, and tensor polarizabilities , were inferred. The brackets 

denote the expectation value of the traceless and symmetric combination of angular momentum generators j 



(^Jih + hji 



We list the results in Table VI and highlight the magnetic moments of and ^He as they coincide within uncertainty 
margins with naive expectations based on the shell model which estimates them as a sum of pp (nn) singlet and 
n {p) moments as recognized in Ref. |82j . This behavior resembles the relations observed at physical 171,^ where the 
deuteron moment is the sum of n and p, while triton (^He) is given by the moment of p (n), approximately. In their 
linear response, the bound nuclei behave similar to external magnetic fields relative to each other. An inference from 
/r = 0 of the a-particle and the moments of the smaller nuclei on the more appropriate cluster model — relative to 
the shell model for the triton and ^He — with ^H-p, ^He-n, and d — d comprising a is: The a ground state resides 
predominantly in the ^He-n and ^H-p configurations with a slightly larger contribution from the former. 

With this data, knowledge about the response of nuclei at large m,r is available. The situation is peculiar enough 
to be rephrased, the measurements considered the coupling of the constituents of nuclei to an external field while 
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disregarding the mutual interaction via the same interaction. In essence, this treatment resembles that of composite 
particles in gravitational fields which takes into account the external field, only, and does not concern itself with the way 
an atom, e.g. , contributes to the bending of space-time. This scenario of large external fields was analyzed in Ref. [H] . 
It was shown that for strong enough fields all two-nucleon bound states at m,r ~ 806 MeV and m,r 450 MeV become 
unbound. At this unbinding field strength the scattering lengths diverge and it was conjectu red^ that at some TOtt this 
field strength is the same for all NN states (see assumption in Ref. [39] introduced in Sec. HI- 

How the spectral data in a small magnetic background field can be utilized by a match to an EFT for the prediction 
of reaction observables was also demonstrated and is summarized below. 

Matching to EFTjf) — Under the assumption that EFT(^) is applicable at ~ 806 MeV, a generalized form 
of the relation Eq. ^ was employecP^to relate energy eigenvalues in the presence of a background field to the LECs 
which couple the magnetic field to nucleons up to next-to-leading order (NLO). The generalization proceeds as in 
Ref. | 86 | and requires the calculations of amplitude poles in a finite volume with the interaction with the gauge field 
at NLO given by (dibaryon version in Ref. [75] , nucleon version in Ref. |7T| ) 

Cihdy = + KiT3 )cr-BN (18) 

ZTTljsf 


and 


/: 2 bdy = [eLi{N^P,Ny{N^P 3 N)B, - P,N)\N^P,N)Bk + h.c.] . (19) 

The electron charge e, the nucleon mass tun, and the nuclear magneton of the neutron (proton) which defines 

the isoscalar (isovector) nucleon magnetic moment Ko(i) = \ (^p(i)^n)j assume values as predicted by LQCD at given 
TO^. The low energy constant Li couples (projector P 3 ) and (projector Pi) and thus contributes at LO to 
the np —>■ dj capture. The operator corresponding to the L 2 LEG does not induce transitions between spin states. 
It contributes to the magnetic dipole moment of the deuteron and is thereby relevanlPTl for asymmetries in the cross 
sections for circularly polarized photons impinging on an unpolarized deuteron target, dy —>■ np. 

In combination with the proper EFT(^) Lagrangian governing the nuclear interaction at NLO, the np —>■ dy 
amplitude can be evaluateci^ . The zeros of the real part of the inverse of this amplitude in a finite volume are related 
to fccot(^^. Scattering lengths and effective ranges for the incoming singlet and outgoing triplet were takeiJ^from 
Ref. |S2| to be degenerate, ^o„p ~ and ri ^ r^. With the isovector ki measured independently in Ref. [52] , Li 
is the only parameter left in the NLO amplitude’s poles which is left undetermined by single-nucleon and scattering 
parameters. It is related to energy shifts between the singlet and triplet eigenstates in the presence of the background 
fielcP^. 

The predictive power was demonstrated by calculating the cross section of the radiative capture with the extracted 
value of Li at TOtt ~ 806 MeVin Ref. [82] and m,r 450 MeV in Ref. [84] where only the value for 806 MeV is quoted: 


(^ 806 inp ^ dy) = 17 (+“ 6 ^) mb . 


( 20 ) 


The asymmetric uncertainty is due to the non-linear input dependence of the cross-section which, in contrast to the 
renormalization of Li, uses different effective-range parameters in the singlet and triplet channel. 

It is noteworthy that the LQCD results of the magnetic moments indicate a similar internal structure of the A = 3 
nuclei as found in nature, i.e., a bound singlet with a single nucleon which determines the spin. Relating this fact — 
the three-body response being given by that of a single nucleon — to the ratio of the separation energy to the binding 
energy of the core singlet, at ^ 806 MeV this ratio is ^ 2.1 compared to ^ 2.9 at physical (see Table VI), the 
structure is expected to change significantly when decreasing the pion mass where the respective ratios at 510 MeV 
and 300 MeV suggest a shallow triton. 


VII. SUMMARY 

The status of a unified description of particle and nuclear physics was presented. This description comprises a chain 
of effective field theories with QCD on the particle, EFT(y-) on the nuclear end, and a bridge through XEFT. 

Interest in the sensitivity of nuclear observables to variations in fundamental parameters — the pion-mass, in 
particular — arose with the inability to solve QCD with the physical pion mass. EFT expansions around = 0 
were employed to assess how sensitive nuclei react on a variation of up to ~ 200 MeV. Those attempts were 
reviewed in Secs. [IT] and [ml The introduced framework of matching a chiral EFT to data at some m^r, and to 
extrapolate nuclear amplitudes to physical or larger m,r while tracing carefully implicit and explicit rrijr dependences. 
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is presented as an ansatz constituting one link in the EFT chain which can be implemented once experimental data 
can be replaced by QCD amplitudes. Work related to the connection between XEFT and EFT(|f), as candidate for a 
few-nucleon EFT, was presented in Sec. [ml 

The following Secs. |IV| through |VI| elaborate on the approximation of QCD in light nuclei by matching contact 
EFTs to lattice calculations. The methodology of LQCD is overviewed with a focus on uncertainties presumably 
responsible for inconsistent EFT postdictions. The available nuclear LQCD data is compiled in Tables 0 and |VI[ 
including only results obtained with controlled approximations. 

In Sec. |V A| and |VB[ we introduced two methods of deriving a nuclear effective interaction. A method based on 
a velocity expansion of a potential consistent with a QCD wave function, and the adaption of an EFT(^) analog to 
LQCD amplitudes. The latter method is included in the summary of efforts to assess the electromagnetic structure 
of nuclei in Sec. EH. 


VIII. OUTLOOK 

It is one aim of this article to review research on the extrapolation of QCD solved at unphysical pion masses to 
physical pion masses as well as the consequences of an enlarged pion mass on nuclear systems. Future work as 
suggested below concerns both aspects. 

a. Outlook in nucleon number — The anecdotal introductory hint to the resemblance between the contemporary 
numerical effort and the historical, experimental one does not apply to larger nuclei. Their properties will not 
be accessible with LQCD within a similar time frame as physical experiments ventured beyond the deuteron. As 
cluster/halo EFT^^ constitute the next link in the EFT chain relating QCD parameters to nuclei with A > 4 and 
have so far been renormalized through a match to datep^ unavailable from LQCD, their connection to EFT(-^) is 
crucial. Whether or not the bulk properties of A > 4 nuclei are a universal consequence of two- and three-nucleon 
data and thus could already be parameterized with m^r is unknown. Specifically, do the unbound Helium-5, the 
shallow aN resonances, the halo structure of Helium-6 emerge at LO EFT(^)? It is an open question if an additional 
renormalization condition in form of a five or six-body counter term are necessary to put those poles in the respective 
LO amplitudes, or if those observables are sensitive to two-body P-wave interactions as which they should not be 
considered before next-to-next-to-leading order. Once EFT and few-body practitioners have addressed this question 
at physical TOtt, the amplitudes can be matched to cluster EFTs and thereby pass the dependence to larger nuclei, 
systematically. 

b. Refining the interaction — The role of the electromagnetic interaction between quarks for lattice nuclei is 
unknown but the response of nuclei to external magnetic fields has been explored in LQCD measurements. The latter 
included a fascinating demonstration how this method can probe extreme conditions inaccessible by experiments. 
These analyses were covered in Sec. |VI| . Next to the numerical effort to implement the electromagnetic interaction 
in LQCD calculations, there remain conceptual issues hamperintheir EFT(y-) consideration. This is an instance 
where EFT can make predictions by including the long-range Coulomb force in EFT(|t) at large pion masses. Contrary 
to nature, the proton-proton system provides both a bound and scattering specimen to study the effect resulting from 
the combination of a long and a presumably relatively short-ranged force. At present, we assume that the energy gap 
between the triton and Helium-3 is approximately an invariant w.r.t. changes in m,r- Noting a peculiar consequence 
of a widening gap, namely a conceivable unbound helion in the presence of a shallowly bound triton, shall motivate 
work in that direction. 

Of interest for understanding the difference between hyper and ordinary nuclei are LQCD measurements at a fixed 
pion mass for both, the SU(3) symmetric point, and with a shifted strange-quark mass rrig, i.e., explicit breaking of the 
flavor symmetry. The effect of, e.g. , an infinitesimal shift in rus on the two states of ^H observed at m,r ~ 806 MeV by 
NPLQCD, could indicate the significance of the SU(3) breaking relative to that of the pion mass for the shallowness 
of the hyper triton w.r.t. to the ordinary triton. The LQCD data available on strange nuclei has not yet been matched 
to a contact theory at 450 or 806 MeV. This requires a generalization of the SU(2) isospin-symmetric Lagrangian of 
EFT(^) to SU(3j^ and thus a comprehensive theory for the baryon octet. In combination with LQCD data, which 
in contrast to nature is roughly as accurate for the s = 0 as it is for the s = — I sector, will allow for a systematic 
study of the peculiar differences between strange and ordinary nuclei. 

The observed m 7 r-insensitivity of the approximate SU(4) symmetry of the nuclear interaction as shown in Fig. 
allows for the investigation of the relevant QCD parameters which cause this remarkable feature. This sensitivity 


general frame work in Ref. m ; treatment of narrow resonances in Refs. |9QI I91| ; applications: proton halJ^, aa in Ref. m . 
two-neutron halo^^, ^Li + n ^ ®Li -h 7 in Ref. HH , ^Be -I- p —> *B -I- 7 in Ref. [96] ^-A^H-ain Ref. EZI ■ 

The work in Ref. ESI can be considered the first rigorous matching of a microscopic and a cluster EFT. 

For recent progress in the systematic treatment of the Coulomb force see Refs. |98tflQQ| . 

We acknowledge the explanation by M. Elyahu and N. Barnea to whom this idea belongs. 
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analysis can be performed at large pion masses since does not seem to be significant for the effect. We are thus 
able to understand the mechanism behind the breaking of Wigner’s SU(4) symmetry also at physical ttItt — an insight 
that would demonstrate the beneficial interplay between LQCD and nuclear EFTs. 

Understanding sensitivity — To predict low-energy renormalization-group fixed points of QCD, we desire 
data about the smooth dependence of nuclear spectra from LQCD. Available lattice calculations investigate nuclei 
at isolated pion masses. For the interpolation between them, no theoretical ansatz is known for m,r beyond the 
convergence rate of ChPT. Even within that radius, the murky formulation of AEFT in the few-nucleon sector 
hampers the uncertainty quantification of those interpolations and renders them less useful. To develop or refine 
interpolating and extrapolating theories, knowledge of whether or not a critical RG flow trajectory is approached will 
be of use. LQCD analyses of energy gaps between nn singlet and the deuteron, deuteron and triton, triton and a, 
and a and Helium-5 at two infinitesimally close TOtt’s would indicate whether or not one approaches a critical for 
any of those systems. 

Such an analysis touches the question how parameters of constituents characterize compounds, e.g. , like the infinite 
two-nucleon scattering lengths in combination with a shallow triton “furnish” the a-particle. The aforementioned 
cluster EFTs work well due to a separation of scales between the excitation energy of the a and the shallow aN poles. 
The ratio between the binding energies of Helium-5 — if bound at all — and the a as one characteristic of nuclei at 
physical m,r, namely their shell structure, is one crucial observable which as a ratio is more accessible to LQCD than 
the bare values. To analyze the ttItt dependence of it would shed light on the emergence of the shell structure, the 
peculiar mass gap at A = 5, and the prominent role of the a as a building block for larger nuclei. 

How nuclear two-body systems can exhibit a peculiar behavior like a Feshbach resonance has been shown by simu¬ 
lating extreme magnetic fields which reside naturally only in cosmological objects like magnetars. Similar features of 
larger systems, like the development of a four-nucleon Efimov spectrum due to a triton near the deuteron threshold — 
a scenario which is admissible within error bars at = 300 MeV — and thus a mass gap at A = 3 is of undeniable 
empirical value to identify the underlying QCD mechanisms for such characteristics of the nuclear chart. To that 
end, ratios of binding energies are more important for the theoretical understanding than relatively less accurately 
measurable absolute binding energies. 

Focal-point system — The five-baryon system as a gateway to heavier nuclei and refined EFTs concludes this ar¬ 
ticle. Firsts it poses a challenge for numerical techniques, for LQCD as well as traditional few-body methods. Second^ 
because its features are neither understood as emergent or unique. Restated, whether or not conventional EFT(|t) ap¬ 
plies to it is unknown. Third, its amplitudes are the canonical candidates for the bridge between single-baryon and 
cluster EFTs. The pion-mass dependence of the dynamics of this system is therefore key to an understanding of the 
emergence of complex phenomena in nuclei from the interactions governing its basic building blocks. 
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